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| ۱ ۱ M" | : ; | surrounds the reactor. The penetration for target injection in the chamber 
1 uud ۲ 5 1 ٍ es | | Wu b | | cns e | | | X | | — م60‎ is very Small and is not considered in this analysis. The total area of 

TE 2 e n cie €. E CEP : , ۱ 5 ; k A ۱ the 20 beam line penetrations in the INPORT tube region is 6 m* implying 

| | 3 ظ‎ that these penetrations occupy only ~ 1% of the blanket volume. This is 

smaller than the statistical uncertainty in the results obtained in the Monte 

: | Carlo calculations used in this analysis. Furthermore, the geometrical ef- 

6* | 8 EL i o T 8 ۱ ۱ Fi 9 ۱ i 3 ۱ | — | fects of these penetrations have negligible effects on the tritium breeding 

ratio and the energy multiplication. Hence, the beam line penetrations are 

not modelled in the present analysis. Detailed radiation shielding analysis 

for the final focusing magnets is given in section V.6. 

: E The neutronics and photonics calculations were performed using the multi- 
Ts rt | m if iij : à ۱ es " di ۱ Ta d i ۱ group three-dimensional Monte Carlo code MORSE.(13) A coupled 25 neutron-21 

gamma group cross section library was used. The library consists of the RSIC 
DLC-41B/VITAMIN-C data 11۳۵ ۰31 and the DLC-60/MACKLIB-IV response data 

library. (4) The combinatorial geometry capability of the MORSE code was used 

to model the problem geometry. Volume detectors were used to estimate the 

here are based on a DT yield of 400 MJ and a repetition rate of 5 Hz yielding SecTIOR AT h 


7.1 x 1029 fusion neutrons per second. 8 point isotropic source was used at 


NO 
Cn 


the center of the reactor cavity with neutron and gamma spectra obtained from 
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3 E du Hn | ۱ ۱ ۱ 2 ۱ ۱ the Monte Carlo calculations yielding less than 2% statistical uncertainties ظ‎ | 6 LL 


in the estimates for the tritium breeding ratio and the energy multiplication. 


| dw ۱ T | | quantities of interest in the different reactor zones. The results presented 
the target neutronics and photonics calculations. 4000 histories were used in y 2 ( 





d ouk Du "f : F Te | | | | | | Ps s 4 Because of symmetry, only 1/60 of the reactor was modeled with reflecting | CROSS SECTIC 


albedo boundaries used at the planes of symmetry. This corresponds to a "pie 


X M —6 008 : f; 0 DET ; | | | TEST | slice" with an azimuthal angle of 6?. The geometry for the computational Figure VI.3-17 HII 


MICROCOPY RESOLUTION TEST CHART | | model used is given in Fig. VI.3-17. To quantify nuclear heating in the 
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The overall tritium breeding ratio is found to be 1.25. The co 
interval for the estimated breeding ratio is 0.025 which is 2X of th 
estimate. The overall tritium breeding ratio obtained here with the 
reactor cavity cylindrical geometry is larger than that obtained fro 
dimensional spherical geometry calculations (1.195). The reason is | 
the spherical geometry case, all source neutrons a: = incident on the 
perpendicularly to its inner surface and hence will see the same blar 
thickness. On the other hand, the source neutrons emitted isotropic. 
the center of the cylindrical cavity iipinge on the inner surface of 
blanket at different angles and hence will see larger blanket thickne 
other words, the neutron source is surrounded by a larger volume of t 
material in the actual cylindrical case as compared to the case of ar 
lent spherical blanket. 

VI.3.3.5 Nuclear Heating 

Table VI.3-12 shows the nuclear energy deposition for neutrons a 
in the different zones. The average power density is also included. 
clear that the contributions from neutron and gamma heating are near] 
same in the breeding blanket while the gamma contribution in the firs 
and reflector is about an order of magnitude higher than the neutron 
bution. This results from gamma generation in the HT-9 structure. A 
of the tota! reactor thermal power comes from gamma heating. The ene 
deposited in the biological shield is 0.06 MeV/fusion which correspon 
power of 6.52 MW. This represents only 0.275 of the total reactor th 
power. Only the energy deposited in the blanket, first wall, and ref 
are considered as recoverable energy. The total recoverable neutron | 
energy in the reactor per DT fusion is 17.553 + .292 MeV which is sli 


less than that obtained for an equivalent spherical reactor (17.95 Mel 
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Table VI.3-11 Tritium Production (Tritons/Fusion) 







Region Zone Number 61 (n,a)T | ^L i (n,n'a)T 


Breeding Blanket 1 0.729 0.018 
0.014 0.00001 


0.212 0.004 


>_> WwW N 


0.235 0.004 


Region Total 1.190 0.026 


Reflector 8 0.022 0.000002 
9 0.001 0.0000001 
10 0.009 0.000002 


11 0.002 0.0000001 
Region Total 0.034 0.000004 
System Total 1.224 0.026 

















dila 


the bottom pool and includes 1 v/o HT-9 structure. The reason in that the 2 
v/o SiC present in the top blanket enhances neutron slowing down and hence 
increases the tritium breeding effectiveness. In fact, our results show that 
the breeding capability of Lij;Pbg4 can be improved considerably by enriching 
Li and/or using moderators in the blanket. The improvement in systems using 
Li17Pb83 15 much more pronounced than that in systems using other breeding 
materials because most of the contribution to breeding in Li; 7Pbg3 case comes 
from the Li (n,a) reaction which has a 1/v cross section in the low energy 


region. 
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Table VI.3-10 Average Scalar Flux (Particles/cm^/Fusion) 


Gammas 


3.237(-8) [.026] 
4.453(-9) [.349] 
2.448(-8) [.056] 
1.108(-8) [.068] 


2.529(-8) [.118] 
1.124(-8) [.600] 
4.328(-8) [.164] 


1.274(-8) [.048] 
3.213(-8) [.245] 
2.230(-8) [.070] 
3.739(-9) [.132] 
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Neutrons 


3.191(-6) [.015]* 
1.691(-6) [.071] 
2.937(-6) [.023] 
1.708(-6) [.043] 


8.302(-7) [.041] 
9,282(-7) [.181] 
1.340(-6) [.047] 


2.384(-7) [.034] 
4.470(-7) [.153] 
3.789(-7) [.044] 
8.354(-8) [.126] 


Region Zone 


Breeding Blanket 1 


First Wall 5 


Reflector 8 





سم 
سم 








“Fractional Standard Deviation 










































VI . 3-6 





vacuum pump, a 2 cm thick region consisting of 50 v/o 316 SS and 50 v/o Cu is 
designated as zone 13 to simulate the cryo-panels. Zones 14 and 15 represent 
inner and outer vacuum regions, respectively. 


VI.3.3.3 Scalar Flux 

Table VI.3-10 gives the average neutron and gamma scalar fluxes in the 
different zones. The fractional standard deviation based on a 68% confidence 
interval is also included. It is clear from the results that the fractional 
Standard deviation is very small in the breeding blanket zones and is rela- 
tively large in the optically thin regions such as the first wall. The re- 
sults show that while the ema flux is about two orders of magnitude less 
than the neutron flux in the breeding blanket zones, it is only one order of 
magnitude less than the neutron flux in the reflector zones. This results 
from the large gamma production following neutron capture in the HT-9 
structure. 
VI.3.3.4 Tritium Production 

Table VI.3-11 shows the results for tritium production per DT fusion re- 
action in the different reactor zones. The contributions from Li (nya) and 
7Li(n,n'a) reactions are shown separately. It is clear that the contribution 
from /Li represents ~ 2% of the total tritium production in the breeding 
blanket region. The reason is that the number of Pb(n,2n) reactions occurring 
in the Li, 7Pbg3 blanket is much larger than that of /Li(n,n'a) reactions。 The 
resulting (n,2n) neutrons are well below the threshold energy for 7Li(n,n'a) 
reaction and hence can produce tritium only through the 61 i(n,a) reaction. In 
the reflector region the Li contribution to tritium production is very small 
because of the neutron spectrum softening in the breeding blanket. 


We notice that as much tritium production occurs in the top blanket as in 


the bottom Li 7Pbg3 pool, even though the top blanket 1s only half as thick as 
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3 ZONE DESCRIPTION 
| INPORT TUBES 
0 BLANKET REGION AT 
VACUUM DUCT 
E 3 TOP BLANKET 
4 BOTTOM LiPb POOL 
® 5 SIDE WALL 
— 3 WALL AT DUCT 
7 TOP LINER 
2 SIDE REFLECTOR 
a REFLECTOR AT DUCT 
0 TOP REFLECTOR 
T BOTTOM REFLECTOR 
2 BIOLOGICAL SHIELD 
I3 VACUUM PUMP 
CROSS 
.ANE 20لا‎ 
=a ۱ ۱ i j 
012345 (m) 
8 REFLECTING BOUNDARIES 





JALL geometry for Monte Carlo calculations. 


== MB حت‎ 
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nfidence | 
ex in ۰ H 
e obtained Table ۷1۰3-12 Nuclear Neat ing reason is that more tritium production is obtained in the cylindrical case | 
Figure VI.3-18 in the biological concrete shield represents only 0.27% of the total reactor 
actual 1 with less neutrons being captured in the HT-9 structure. A neutron absorbed 1 
, thermal power. The thermal م‎ for the HIBALL power plant is 10,193 MW(th). 
n the one- ۱ — —— ——— — in 9Li releases ~ 4.8 MeV while if it is captured in the HT-9 structure ~ 7 | | ۳ "— p p 3 M( th) 
1 a, MeV /fusion Density (W om VI.3.4 Time-Dependent Neutronics and Photonics 
hat in — P"— MeV is released. | ENERGY FLOW IN HIBALL 
i im d | HIB PULSE ; VI.3.4.1 Introduction 1 
blanket Neutrons Ganma s In order to get a statistically adequate estimate of nuclear heating in (4.8 MJ) 
١ ۱ ۱ ١ In an inertial confinement fusion reactor the neutron source has a pulsed 
ket the vacuum pump with a reasonable number of histories, a Monte Carlo calcu- 
ing 1 46 49 4.409 i i — nature because of the very short burn time over which the fusion reactions 
lly at Breeding Blanket lation with an angular source biasing was performed. In this run, 4000 d 
2 0.074 0.017 3.911 ۱ R -11 3 TARGET occur (10-100 ps). Furthermore, the neutron pulse does not reach the first 
the histories were used. A volumetric nuclear heating rate of 6 x 10 W/cm” was 17.6 MeV/FUSION 
A 0 e haiti i surface of the blanket until ~ 100 ns after the burn and the neutron slowin 
5565۰ In " — e — obtained in the vacuum pump (zone 13). The fractional standard deviation was ۱ 400 MJ/PULSE 9 
0 0 e 1 i ; rá | down time in the blanket is much greater than the duration of the neutron 
reeding ۰ 9 — — 0.25. Because the vacuum ducts do not see direct line-of-sight source 9 
i ; 0 ° 1 source. Significant softening of the energy spectrum of neutrons escapin 
| equiva- — er ا‎ — X ( neutrons and they are bent twice, neutron streaming through the ducts was 9 9 gy SP ping 
| from the target results from neutron interaction with the dense (~ 102°/cm?) 
First Wall 5 0.004 0.050 1.222 found to not cause any serious problem to the vacuum pump. 
| R — à; fuel and surrounding tamper material. This leads to a considerable time of 
| 6 0.0001 0.0004 2.010 The energy flow for the HIBALL fusion reactor design is illustrated in NEUTRON & GAMMA X-RAY & DEBRIS LOST IN ENDOERGIC REACTIONS 
۱ ۱ ۰ ۱ 12 559 MeV flight spread as neutrons reach the first surface. Therefore, time dependent 
nd gammas ۱ 7 0.003 0.039 3.068 Fig. VI.3-18. The values given for the power correspond to one reactor cavi- . e 4.862 MeV 0.179 MeV | | 
NEAN Restate th 285.42 MJ 110.49 MJ 4.08 MJ | neutronics studies are essential for the proper analysis of inertial confine- ۱ 
It is Region Total 0.007 0.089 1.653 ty. Therefore, these values need to be multiplied by 4 to calculate the power 1427.1 MW 552 MW 20.35 MW Sul Manes نی‎ ۱ 
y the from the whole power plant。 This corresponds to a total power plant thermal à 
° 1.628 0.939 As a result of the pulsed nature of the neutron source, high instantane- 
t wal! Reflector ۲ — power of 10,193 MW(th). The overall energy multiplication defined as the / 8 
9 0.007 0.066 1.020 i ١ ۱ ous damage rates are present in an inertial confinement fusion reactor wall 
contri- total energy deposited in the system, including the energy deposited by X-rays | NEUTRON & GAMMA X-RAY & DEBRIS TOTAL 
10 0.070 0.733 1.465 - ar BREEDING 14.659 M and structure. This can lead to significant changes in the microstructure (14) 
bout 60% and target debris at the first surface of the blanket, divided by the fusion eV 4.862 MeV 19.521 MeV 
11 0.013 0.126 0.257 | BLANKET 333.16 MJ 110.49 MJ 443.65 MJ of the first wall material. It has also been found that high instantaneous dpa 
rgy reaction yield of 17.6 MeV, is found to be 1.274. 1667.5 MW 559 MW 2219.5 MW (87.05%) 
Region Total 0.245 2.553 0.941 | : ` rates result in higher recombination rates with the void growth being in- 
ds to a VI.3.3.6 Summary ie ns) 
hibited and swelling decreased. Accurate instantaneous damage rates can 
ermal System Total 7.901 9.652 3.345 A three-dimensional Monte Carlo neutronics and photonics analysis was 0.096 MeV 0.096 MeV 22.415 MeV 3 9 
| 8 ۱ ۱ WALL 2.24 MJ 2.24 MJ 509.39 MJ be calculated by performing time-dependent neutronics studies. Previous time 
lector performed for the HIBALL reactor cavity. An overall tritium breeding ratio of 10.8 MW 10.8 MW (0.42%) 2548.3 MW 
er 559 dependent neutronics studies have been used to calculate the instantaneous dpa 
and gamma 1.25 and an overall energy multiplication of 1.274 were obtained. The tritium 
2.798 MeV 2.798 MeV rate in the first wall of a laser fusion reactor(16) and an electron beam 
ght ly breeding ratio is higher than that for the equivalent spherical reactor cavi- REFLECTOR 63.5 MJ 63.5 MJ 
ERAM : : fusion reactor. (17) No correction for the time of flight spread of neutrons 
/). The ty. Nuclear heating in the vacuum pump was found to be very small. The power 318.0 MW 318.0 MW (12.53%) 
within each energy group was made in these studies. 







OVERALL ENERGY MULTIPLICATION = 1.274 
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dependent blanket neutronics analysis performed using the time dependent 
multi-group discrete ordinates code TDA. ANISN was used to calculate the time 
integrated reaction rates. Since the particle transport codes used are one- 
dimensional, spherical geometry was used in the blanket calculations and hence 
the results represent the worst conditions at the central plane of the 
cylindrical reactor. A P3-S8 approximation was used in the transport calcu- 
lations. A coupled 25 neutron-21 gamma group cross section library has been 
used. This library consists of the RSIC DLC-41B/VITAMIN-C data library(3) and 
the DLC-60/MACKLIB-IV response data library. (4) 

VI.3.4.3 Neutron Source 

The target, though absolutely very small, is an extremely dense medium 
composed primarily of light elements. This results in a substantial collision 
probability for fusion neutrons created within it. The 14.1 MeV fusion 
neutrons are degraded in energy as a result of elastic collisions with the 
fuel (DT) and inelastic collisions with the high Z materials (Pb). Neutron 
multiplication also results from (n,2n) and (n,3n) reactions with the con- 
stituent elements of the target. Since the neutrons escaping from the target 
represent the neutron source for the blanket neutronics calculations, it is 
essential to perform detailed neutronics calculations for the target to 
account accurately for spectrum softening and neutron multiolication. 

Since the burn time of the target (10-100 ps) is much greater than the 
slowing down time of fusion neutrons in the extremely dense target )- 0.1 ps), 
steady state calculations yield useful information. Furthermore, since the 
burn time is much smaller than the time neutrons take to reach the first sur- 
face of the blanket including the time of flight spread, neutrons are assumed 
to be emitted in a pulse with zero time duration (6(t)). This implies that 


the time integrated spectrum of the neutrons leaking from the target as calcu- 








A modified version of the time dependent discrete ordinates code TDA(18 ) 
has been used to perform time-dependent neutronics analysis for HIBALL. The 
effects of the INPORT first wall protection concept on the peak instantaneous 
and average dpa and gas production rates in the ferritic steel first wall have 
been investigated. The Steady state discrete ordinates ANISN (1) code was used 
to determine the average time integrated radiation damage rates. 

Energy deposition in the porous tubes and first wal] of the HIBALL fusion 
reactor is important in the determination of their structural integrity. The 
temporal and spatial dependence of energy deposition is also required for 
stress analysis studies. The energy deposition rate integrated over the vol- 
ume of the blanket of a laser fusion reactor was calculated in a previous 
study. (19) However, no detailed time dependent fluxes or energy deposition 
rates at different points were given. In this work, the instantaneous energy 
deposition rates at different positions in the blanket and first wall are 
calculated. 

V1.3.4.2 Heavy Ion Beam Reactor Computational Model 

A schematic of the blanket, first wall, reflector, and shield configu- 
ration for the HIBALL fusion reactor is given in Fig. VI.3-1. The nuclide 
densities used in the calculations are given in Table VI.3-1 The results 
presented here for the damage and energy deposition rates are based on a DT 
yield of 400 MJ which corresponds to 1.42 x 1020 source neutrons per pulse. A 
repetition rate of 5 Hz is used to determine the total dpa and helium produc- 
tion in the first wall per full power year (FPY). 

The steady state discrete ordinates code ANISN was used to perform de- 
-ailed neutronics and photonics calculations within the target, giving the 
ime integrated energy spectrum of neutrons escaping from the target as ex- 


lained in section 111.1. This Spectrum was used as a source for the time 
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Figure VI.3-20 Arrival 
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As neutrons travel from the target to the first surface of the blanket 
considerable time of flight spreading occurs because of the broad energy 
distribution of these neutrons. The arrival time spectrum at the first sur- 
face located 5 m from the neutron source is shown in Fig. VI.3-20. The num- 
bers at the bottom of the graph indicate the energy groups associated with the 
various times. Since all neutrons in a specified energy group are forced by 
the multigroup treatment of the neutronics code to travel at the same velocity 
corresponding to the average group energy, these neutrons will arrive at the 
first surface simultaneously and no time of flight spread corresponding to the 
energy group width will be observed. In order to preserve the correct arrival 
time spectrum at the first surface, source neutrons are given artificial birth 
times associated with their energies within the group. If the time of flight 
spread corresponding to the gth energy group width is Ag， the modified neutron 
source for group 0 is represented by a rectangular pulse with a time duration 
of Ag. This necessary correction for time of flight spread of neutrons within 
each energy group was not made in the previous time dependent neutronics 
analysis of inertial confinement fusion reactors. 

VI.3.4.4 Time Dependent Neutron Spectrum in the HT-9 First Wall 

A version of the time dependent multigroup discrete ordinates code TDA, 
modified to facilitate its use for the analysis of inertial confinement fusion 
reactors (20), has been used to perform time dependent neutronics and photonics 
calculations for the HIBALL blanket and reflector model illustrated in Fig. 
VI.3-1. For comparison, calculations have been also made for the case of an 
unprotected wall. 

The source for the problem involves Dirac delta functions which are quite 
difficult to represent with standard finite difference methods. To circumvent 


this difficulty the uncollided flux is determined analytically and an analytic 
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lated by the steady state code ANISN gives the correct neutron source spectrum. 


In other words, if n (E,t) is the exact spectrum of neutrons leaking from the 





target at time t within the burn, the time integrated spectrum calculated by 
























ANISN is given by IO ۲ 3 FY TOTO VW NN FNN 2 ' : 
- SPECTRUM OF NEUTRONS N 
Latii: s xu as — - LEAKING FROM HIB TARGET || - 
Since the burn time is very small, all neutrons a 1 


where tp is the burn time. 


leaking at different times during the burn are assumed to be emitted at t=0 










and the source for the blanket calculations is given by 
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S(E,t) = n, (E) 6(t) . (VI.3-2) 
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| 





Since the target has a very small radius compared to the radius of the reactor 





cavity, the neutron source is represented by a point source at the center of 


< 
N 


the reactor. Furthermore, neutrons are emitted isotropically from the target 





NEUTRONS LEAKING/FUSION 





and the complete representation of the neutron source is given by 









S(u,r,E,t) = n, (E) Slr) 6(t) (w1) . (VI.3-3) 
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Detailed neutronics and photonics calculations have been performed for 









the HIBALL target as discussed in Section III.1. The results show that a 










target neutron multiplication of 1.046 is obtained. This results mainly from 







(n,2n) reactions in the dense DT fuel core and LiPb pusher. The average ener- 






Fig. VI.3-19 Spectrum of neutrons leaking from the pellet. 


(The ordinate is the integral in each 
energy group.) 


gy of neutrons leaking from the target is 11.98 MeV implying that 71.2% of the 






fusion energy is carried by the emerging neutrons. The spectrum of neutrons 






leaking from the target, n (E), is given in Fig. VI.3-19. 


VI.3-61 




















BM سد‎ 


with a very small energy loss per collision resulting in much longer sl 
down times. Therefore, the time spread of the spectrum in the first we 
our system is determined primarily by the slowing down time in the blar 

On the other hand, the time of flight spread is the dominant fact: 
the case of the unprotected wall because of the absence of the "inner" 
down region. The spectrum at 5 ns after the arrival of the leading ed 
the neutron pulse consists mainly of first group neutrons coming direc! 
the source. The soft part of the spectrum shown in Fig. VI.3-22 is du 
backward scattering of 14.1 MeV neutrons in the reflector region and ne 
slowing down in the wall itself. As time elapses, lower energy source 
neutrons arrive at the first wall. The effect of neutron slowing down 
reflector is not as pronounced here as compared to the effect of neutr 
slowing down in the inner blanket for the case of protected wall. The 
is that the probability of backscattering of high energy neutrons is S 
compared to the probability of forward scattering. Therefore, the tim 
dependent damage rate in the unprotected wall is affected primarily by 
time of flight spread. 
VI.3.4.5 Atomic "isplacement Rate 

The instantaneous damage rates in the first wall have been calcul 
using the time dependent neutron spectrum in the wall and the appropri 
action cross sections. The instantaneous dpa rate in the protected fe 
steel first wall is given in Fig. VI.3-23. The results correspond to 
fusion yield. The cumulative dpa is illustrated in Fig. VI.3-24. It 
that the damage occurs over a relatively long time resulting in a peak 
taneous dpa rate of 0.009 dpa/s at 140 ns after burn. This peak corre 
to the 14 MeV source neutrons arriving without collision. A broad pea 


~ 270 ns corresponds to neutrons emitted in (n,2n) reactions with lead 
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Figure VI.3-22 Neutron spectrum in the unprotected first wall. 
(The ordinate is the integral in each energy 


group. ) 
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Figure VI.3-21 Neutron spectrum in HIBALL first wall. 


(The ordinate is the integral in each energy 
group. ) 


0 


Io? 





NEUTRONS/cm? /s / FUSION 








VI.3-50 — 


first collision source is used by the code to determine the collided part of 
the flux. 

The neutron spectrum in the protected ferritic steel first wall is illus- 
trated in Fig. VI.3-21. The results are shown at 10 ns and 70 ns after the 
uncollided neutrons of the highest energy strike the wall. The neutron 
spectra in the unprotected wall at different times following the arrival of 
the leading edge of the pulse at the wall are given in Fig. VI.3-22. Com 


paring the results in Figs. VI.3-21 and VI.3-22, we notice that significant 


neutron attenuation occurs in the INPORT tubes. 


Neutrons slowed down in the inner blanket region take a relatively long 
time to reach the first wall. This time depends on the slowing down time of 
neutrons in the blanket, the speed of slowed down neutrons, and the position 
within the blanket at which the slowing down interaction occurs. At early 
times, following the arrival of the fastest uncollided neutrons at the first 
wall, most of the neutrons in the first wall are in the high energy groups. 
Only neutrons which have been slowed down in the porous tubes close to the 
first wall are able to arrive at the wall by this time. Neutrons slowed down 
in the wall itself result in the relatively low peak at - 1 MeV in Fig. 
VI.3-21. At 60 ns later, most of the high energy neutrons have already passed 
the first wall while neutrons slowed down in the blanket are still arriving at 
the first wall primarily because of the relatively long slowing down time of 
neutrons in the blanket. Therefore, the time dependent spectrum and conse- 
quently the instantaneous damage rates in the first wall are governed by two 
factors; the time of flight spread and the slowing down time spread in the 
blanket. In lead, which is the main constituent of the blanket, the slowing 
down time from 14.1 MeV to the inelastic threshold energy of 0.57 MeV is about 


l us. At energies below 0.57 MeV, slowing down is due to elastic scattering 
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Figure 1.3-25لا‎ DPA rate in the unprotected first wall. 
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INPORT tubes. This considerable time spread results from the relatively long 
slowing down time in the blanket allowing neutrons of energies greater than 


the dpa threshold energy of iron )- 1 keV) to exist in the first wall over a 


long period of time. 
The instantaneous dpa rate in an unprotected wall is shown in Fig. VI.3-25. 


The time spread is very small compared to the case of the protected wall. The 
time spread is determined here by the time of flicht spread of neutrons as 
they travel from the source to the wall. The effect of the INPORT tube pro- 
tection, used in HIBALL, on the atomic displacement in the first wall is given 
in Table VI.3-13. The total dpa per full power year (FPY) was determined 
using the steady state code ANISN with a repetition rate of 5 Hz. 


The wall protection is found to decrease the total cumulative dpa and the 
peak instantaneous dpa rate by factors of 9.4 and 1190, respectively. The 
larger reduction in the peak instantaneous dpa rate results from time spread 
FWHM increasing from 5 to 980 ns. The reduction in dpa achieved here is much 
larger than that achieved in a 316 SS first wall protected by 0.5 m of liquid 
Li 07) typical of the HYLIFE concept. (21) In this case the average rate was 
found to decrease by a factor of 3.6 and the peak instantaneous rate to de- 
crease by a factor of 13. Larger reduction is obtained in our case because Pb 
is more effective than Li in slowing down high energy neutrons and C is more 
effective than Li in slowing down low energy neutrons. If we assume that the 
life of the wall will approach 0 Mw-y/m* (the target for the U.S. fusion ma- 
terials program), and 1 MW-y/m = 11 dpa, then only 4.9 MW-y /m- will be ac- 
cumulated with the INPORT units and 46 MW- y/m2 without them. This implies 
that the wall protection concept used in HIBALL will be very effective in re- 


ducing damage to the first wall, allowing it to last for the lifetime of the 


reactor (20 full power years). 
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Figure VI.3-24 Cumuiative DPA in the protected first wall. 
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Figure VI.3-23 DPA rate in the protected first wall. 
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Figure VI.3-26 Helium production rate in the protected first 
wall. 
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Table VI.3-13 
Effect of First Wall Protection on Atomic Displacement 


Peak Instantaneous Time Spread fwhm Total dpa 


dpa rate (dpa/s) (ns) (dpa/FPY) 
Protected Wall 0.C09 980 2.7 
Unprotected Wall 10.7 5 25.4 


VI.3.4.6 Helium Production Rate 


The instantaneous helium production rate in the protected ferritic steel 
first wall is given in Fig. VI.3-26. A peak instantaneous helium production 
rate of 0.11 appm/s occurs about 15 ns after the leading edge of the pulse 
arrives at the wall. The cumulative helium production is given in Fig. VI.3-27. 
t is clear that the time spread here is much smaller than that for the dpa. 
he reason is that the (n,a) reaction in iron has a threshold energy of ~ 2.7 
leV and neutrons at energies greater than this energy exist in the first wall 
ver a relatively short period. Only one peak occurs because the (n,2n) 
eutrons do not contribute to helium production. The instantaneous helium 
roduction rate in an unprotected wall is given in Fig. VI.3-28. Again the 
ime spread in this case is determined primarily by the time of flight spread 
f source neutrons. 

The effect of the INPORT first wall protection on the helium production 
5 given in Table VI.3-14. The total helium production in a full power year 
5 found to decrease by a factor of 630 whi le the peak instantaneous helium 


roduction rate is found to decrease by a larger factor of 1627. Since the 


lium production reaction is a high energy reaction, neutron slowing down in 
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Table VI.3-14 


Effect of First Wall Protection on Helium Production 


Peak Instantaneous Total Helium 
Helium Production Time Spread FWHM Production 
Rate (appm/s) (ns) (appm/FPY) 

Protected Wall 0.11 26 0.364 

Unprotected Wall 179 5 737 


the blanket is found to have a pronounced effect on helium production in the 







wall. 


VI.3.4.7 Energy Deposition Rate 





The time-dependent neutron and gamma fluxes are used together with the 


appropriate kerma factors for neutron and gamma energy deposition to calculate 





the energy deposition rate in the blanket and first wall. The results are 






given in Fig. VI.3-29 at the blanket first surface, the center of the blanket, 








and the first wall. The time distribution is very narrow at the first surface 





and broadens as one moves towards the first wall. While the time spread at 


















the first surface is determined by the time of flight spread, the spread at 
the first wall is determined by the slowing down time in the inner blanket. 













For a 400 MJ fusion yield the peak instantaneous power densities in the 
blanket and the first wall are found to be 1.82 x 108 and 2.73 x 10° W/cm, 





respectively. This corresponds to peak to average temporal power density 






ratios of 8.48 x 106 and 1.65 x 10?, respectively. The results in Fig. 






VI.3-30 give the instantaneous energy deposition rate in the unprotected first 






wall. The distribution is very narrow with a peak power density of 3.3 x 108 





Helium production rate in the unprotected first wall. 
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Figure VI.3-28 
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Figure VI.3-27 Cumulative helium production in the 
protected first wall. 
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VI.3-63 VI.3-64 VI.3-65 | 
3 71.3.5 Radioactivity and Afterheat FIRST 
W/cm" for a 400 MJ fusion yield. The peak to average temporal power densit NKET ALL | REFLEC. SHIELD 
— -— Radioactivity will be induced in the coolant and structure of the HIBALL — — - ۳ 
ratio in the wall is found to be 1.1 x 107. The INPORT concept is found to ملد‎ E e 
reactor through activation by the fusion neutrons. To treat the probl em 
Nh, ee ae www ا ل‎ decrease the peak instantaneous power density in the wall by a factor of RVALS 1 1۲. 20 INT. 10 INTERVALS 
3 exact1y energy dependent neutron fluxes must be calculated for each point in 0 INTERVALS 5 INTE 
d ~ 1210 and the total nuclear heating in the wall by a factor of ~ 18. These ١ ۱ 57 
ren T T ۱ T — the cavity. Since the ability to perform these calculations with reasonable * 500.0 7000 701.0 741.0 
7 1 results are useful for stress analysis studies. | 
a= as an am —* BLANKET E pe ع‎ E ME 3 cost and effort does not exist at this time a simplified approach was taken. cm 
۰ L VI.3.4.8 Summary and Concluding Remarks ١ 
— — ۲۶7۲۱ FIRST WALL 一 1 1 تنس‎ Specifically the reactor was modeled in spherical geometry with materials and | COMPOSITIONS | 
3 R = Tm | Time dependent neutronics analysis for the ferritic steel first wall of 
7 " ۱ ۱ dimension consistent with a cut through the midplane of the reactor. The 
i lo? ۲ a the HIBALL fusion reactor design has been performed. The analysis accounts a ۱ ۱ 
d ' | multigroup neutron flux used in the activation calculations was taken from the BLANKET FIRST WALL . REFLECTOR SHIELD 
1 : 4 for neutron interactions in the target resulting in neutron multiplication and 1 i ١ ۱ i 
: d one-dimensional ANISN calculations described in section VI.3.1。 For the pur- 
: < significant spectrum softening. The time dependence of the neutron source is ۰ 4 HT-9 13 w/o Fe 
一 = | d | | pose of the calculation the reactor was broken up into 37 intervals and the ۱ 99.3 w/o Li17Pbg3 100% HT-9 88 * » —“ 
4 | eo modified in such a way that the multigroup treatment adopted in the time de- | ۱ | Í ۳۲ | 0.7 w/o SiC 12 w/o 17 83 
Mi 1 0 ~ average flux in each interval was used. The size and compositions of the .5 w/o H 
5 » 20 pendent transport code predicts the correct time of flight spread of neutrons n ۱ é "appe (22) 44 w/o O 
N - u d | various regions is shown in Fig. VI.3-31. The DKR code was used to calcu- 7 w/o Ca 
5 P in each group as they travel from the source to the first surface of the i 96 TM ۱ ۰ i 
~ | ~ late the radioactivity parameters -- activity, afterheat, and biological 4 w/o Al 
— e blanket. A modified version of the time dependent discrete ordinates code TDA DENSITIES 3.5 w/o Oth 
d ۰ n 4 ~ 106 一 5 hazard potential (BHP). 
~w نت جه‎ ee w) d has been used. 3 
۱ 3 E - The activities were calculated for an operating time of two years. Pre- LiPb: 9.4 g/cm 
NS 4 2 E Neutron slowing down in the INPORT first wall protection system is found ١ ۱ ۱ ۱ ۱ SiC: 3.17 g/cm? 
4 — d `~ | vious calculations on other fusion reactor designs show that the buildup of 3 
À — © > to have a significant effect on the time dependent spectrum and damage in the f PI i ۱ ۱ PTEN HT-9: 7.8 g/cm 
EJ radioactivity is quite rapid, becoming significant after only à few hours of | 
: 5 first wall. The time over which the damage occurs is found to be determined ۱ : . | Figure VI.3-31 
3 operation. By the end of two years all short-lived 1sotopes saturate and the 
0 ١ primarily by the slowing down time in the blanket. In the case of an unpro- | 
| ~ ۲ P buildup continues at a slower rate. Thus the activity at times after shutdown 
< tected wall, where no slowing down occurs in front of the wall, the spread is 
. 5 105 ۱ ۱ 1 ۱ ۱ ' - of much less than two years is relatively independent of operating time while 
«i determined primarily by the time of flight spread. 
| |20 130 140 150 160 I7O 180 — the activity at times very much longer than two years after shutdown is ap- 
a a ee A mo s 4 4 Using the INPORT first wall protection concept results in significant re- 
i30 50 170 190 210 230 TIME FROM BURN (ns) / " proximately proportional to operating time. 
ductions in peak instantaneous and total dpa and helium production rates al- —-— i i i 
The total activity per unit of thermal power for the cavity region, first 
FROM BURN (ns) ۱ Figure ۷۲۰3-30 Energy deposition rate in the unprotected first lowing the first wall to last for the reactor lifetime (~ 20 FPY). Our re- j ۱ —— 
wall. material wall and reflector is shown in Fig. VI.3-32. The level at shutdown 






sults also show that the peak power density in the first wal! resulting from 





deposition rate in HIBALL blanket anc first is 0.62 Ci/watt. this value is somewhat larger than might be expected con- 








nuclear heating decreases considerably when the porous tube concept is used to ۲ ۱ E ۱ ۱ ۲ 
sidering that the bulk of the activity is due to the steel which is shielded 






protect the wall. 
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Figure VI.3-34 


Figure VI.3-33 
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from the neutron source by the Pb-Li INPORT region. However, the neutron 
multiplying effects of the Pb in the INPORT region results in the flux in the 
reflector region remaining relatively high and this feature combined with the 
large amount of steel in the reflector produces the high activity level. As 
is typical with steel systems the activity falls off rather slowly with time 
after shutdown requiring approximately three weeks to be reduced by a factor 
of 10 and two years to be reduced by a factor of 100. The relative contri- 
butions of individual isotopes is shown in Fig. VI.3-33. The initial activity 
is dominated by the contributions of 25mm and 203pb。 56Mn (T1/2 = 2.58 hr) 
decays away first followed by “03pb (T1/2 = 52 hr). Note that unless flow 
were maintained, the 203 pb activity would not be present in the reactor 
itself but would be associated with the coolant storage facility. The next 
major decrease in activity is due to the decay of م51‎ (T1/2 = 27.7 days). It 
is followed by the decay of SFe (T1/2 = 2.7 yr) which produces the large 
activity drop between 1 and 100 years. The long term activity is due to 9 3Mo 
and 93Nb originating in the small amount of Mo in the steel (HT-9). 

The afterheat is shown in Fig. VI.3-34. The value at shutdown is 0.667 
of the operating power. For the same reasons as the radioactivity this value 
is somewhat higher than expected, but is consistent with the activity level. 
Since the decay energies of the various isotopes are different the afterheat 
curve has a shape which differs from that of the activity. It falls off much 
more rapidly, by a factor of 10 in about two days and a factor of 100 in two 
weeks. The re'acive contributions of the various isotopes are shown in Fig. 
VI.3-35. At shutdown 56Mn is the dominating isotope. The decay of this iso- 
tope plus the relatively rapid decay of 203pb account for the order of magni- 
tude drop in two days. ler, which contributes significantly to the activity 


in the one week to one month period, contributes very little to the afterheat, 
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thus the afterheat in this period is governed primarily by the decay of 203pp. 


The decay energy of Fe is also low as is the decay energy of the remaining 





isotopes which leads to the rather low residual afterheat subsequent to about 















































r ۱ ' | 1 + ， ۲ r 100 — ۱ ۱ ۱ 
IO ۲ FRACTION OF THE TOTAL AFTERHEAT - The Biological Hazard Potential (BHP) is shown in Fig. VI.3-36. BHP is 
OF THE MAJOR CONTRIBUTING ISOTOPES defined as the ratio of the activity present in the system per unit of power 
0.9 F- HIGALL “203 y^ giren Mo93 ١ to the level of activity allowed per unit volume of air in the U.S. Nuclear 
56 Regulatory Commission Regulations (10 CFR 20). Thus it is a measure of the 
0.8 F- = potential hazard of a radioactive material. The BHP of the coolant first wall 
d and reflector is shown in Fig. VI.3-36. The units are cubic kilometers of air 
d per kilowatt of ther»al power. The shape of the curve is similar to the shape 
ce L- 一 | of the activity curve being 26.0 at shutdown and requiring almost one month to 
be down by a factor of ten and with a long term reduction of only three orders 
0.5 F- P of magnitude. ۱ 
The activity, afterheat and BHP of the shield are Shown in Figs. VI.3-37, 
0.4 [- 1 VI.3-38, and VI.3-39. The neutronics calculations for the shield were per- 
4 | formed for a concrete with no reinforcing material. For the purpose of ob- 
oor taining a better estimate of the radioactivity related parameters of the 
0.2 上 shield the fluxes from the above calculations were used with a concrete having 
L. — reinforcing naterial and with composition typical of the sacrificial shield of 
O. F im 1Om Ih 6h lid Iwk Imo ! y ۱0۴ ۱00 000 1 a BWR. (23) It is this composition which was shown in Fig. VI.3-31. The 
۱ a A | 1 ۱ T activity in the shield is significantly lewer than the activity in the reflec- 
0 io! 102 io? 4 5م‎ 05 IO7 io8 10? IO id! tor, but the value of 6.3 x 1073 Ci/watt is still significant. The time be- 


havior of the activity is dominated after the first day by the activity of 


TIME AFTER SHUTDOWN (sec) 


55Fe which continues to be the major contributor until times greater than 100 


Figure VI.3-35 years. The afterheat at shutdown is 7.7 x 10-29 of the operating power which 


while relatively small, will stil] require some residual heat removal capacity. 
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Figure VI.3-38 
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Figure VI.3-37 
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beam transport system. The omission of these details do not invalidate the 
calculations presented but do require that they be used with care. 

The calculations are also only as good as the cross section data used in 
them. While the data set used in DKR(24) was compiled with care using the 
best available data either experimental or calculational, it is likely that 
some reactions have been omitted. For example it is known that if lead is 
exposed to a neutron flux, an isomeric state at 0.91 MeV in 204pb is excited 
(T1 /2 - 66.9 min). However, the cross section for this reaction is not known 
and therefore has not been included. 

As a final point the effect of impurities has not been accounted for in 


the LiPb. Bismuth, for example, is a common impurity in commercial lead. The 


reaction 2095; + N + 210p; £ 210po is then possible and of concern because 


is alpha emitting and volatile. This problem was considered in the‏ 0م210 
NUWMAK (25) study. A Bismuth concentration of 42 appm was used as a basis for‏ 
the calculaticn. This value was taken from the assay of a relatively high‏ 
purity lead purchased for a slowing down time spectrometer. Assuming Bismuth‏ 
content was calculated to be‏ 6م210 has an 1/v absorption cross section the‏ 
Ci/tonne of Pb. There is approximately 2300 tonnes of Pb/cavity in‏ 0.26 
HIBALL. If the previous ratios were correct the corresponding amount of “10po‏ 
would be ~ 600 Ci. Inis amount of activity, while significant, is very much‏ 
less than the activity of the rest of the material in the system(Should pre-‏ 


sent no special problem. Also, since the half-life of 210pp is short (138.4 


days) it oresents no long term storage problem. 
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The afterheat drops off faster with time than does the activity primarily 
because of the low decay energy of 55Fe。 The BHP is also relatively low and 
falls off slowly, again because of the hazard characteristics of nuclei in- 
volved. To get a somewhat different appreciation of the magnitude of the 
activities in the shield the above values at shutdown correspond to an average 
activity of 3.4 x 1073 61/603 , an afterheat of 4.1 x 107? W/cm, and a BHP of 
6.7 x 10-11 km3/cm3. 

The calculated neutron and gamma doses outside the shield have been pre- 
sented in section VI.3.2. An additional component is the dose due to acti- 
vation products. This component is present during operation and shutdown 
periods. To calculate this dose the activities calculated above have been 
used along with a one-dimensional gamma transport program. The results of 
this calculation indicated that the dose outside the shield from activation of 
the reflector and materials in the shield is quite low, 4.7 x 107? mrem/hr at 
shutdown. 

These low values plus the consideration that the neutronics and prompt 
gamma calculations were done with no high Z materials in the shield, reinforce 
the idea that the shield is somewhat overdesigned and that a more realistic 
choice of shield composition would lead to a thinner shield using less materi- 
als while still showing acceptable dose levels outside the chamber. 

These calculations must be used with a certain amount of care. First the 


geometric modeling of the system is very simple. It would have been much more 


satisfactory to have three-dimensional fluxes for the actual geometry. Further- 


more no account has been taken of some of the special features of the system. 
For example, activation and shielding of the vacuum system has been ignored. 
The activation of the tube support structure in the INPORT region has been 


omitted. Aiso no calculations have been performed to get activities in the 
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Figure VI.3-39 
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VI.4-1 
VI.4-2 — VI.4-3 
n Center, VI.4 Behavior of Pb-Li Vapor | | 
A unique engineering problem for ICF reactors is the design of a first vaporization of the Pb-Li is described. Section VI.4.2 deals with the where & 5 energy/unit volume/unit frequency © 
surface which can carry away the steady-state average surface heat and at the stopping of target generated ions in the vapor, the flow of Pb-Li into the Foy = incident intensity, energy/unit area/unit frequency 
rt ۰ : 
1 same time survive the transient temperature increase generated by X-ray and reactor chamber and the radiation of photon energy from the vapor back onto u, = absorption cross section at energy hv. o 
ion debris from the target explosion. In HIBALL, the INPORT tubes serve this the tubes. Section VI.4.3 contains a description of the vapori zation and The total energy deposited per unit volume is obtained by integration < 
: ۱ wap Qc 
purpose. The INPORT tubes, filled with flowing Pbg3li,7 coolant, protect the condensation of the Pb-Li vapor. The pressure history of the cavity 5 thus over the X-ray spectrum. INTEGRATED RADIATION SPECTRUM | : 
Code for first structural wall from excessive neutron damage. The SiC tubes themselves established. Different X-ray spectrums result in different masses of Figure VI.4-1 shows the X-ray energy spectrum obtained from the target 26 
170, ۱ ۱ ١ ۱ ١ ۱ ۱ 
are protected from the short range X-ray and ion debris by a thin layer of evaporated material and this is also investigated. calculations (3) discussed in section 111.1. The spectrum 1s peaked around 3 Time = 3.5109 Nanosec B 
logy» PbLi that flows down the outside of the porous tube. VI.4.1 X-ray Energy Deposition and the Resulting Evaporation keV and extends up to 300 keV. Since most of the photons have energies well lO p"rr9£9.v v9." "ad FEWw CU. 
ctor E , ۱ 
rthwest After a target explosion, the X-ray energy is deposited within 10-3 cm of Target generated X-rays will deposit their energy over 3 short range in below 100 keV, we have used the Biggs data. s o 
0. ۳۳ ۱ b ۹ ^ 
* | the first surface. The temperature of a thin layer of Li17Pba4 exceeds its Pb-Li and will heat this thin layer to high temperatures. The X-ray depo- Figure VI.4-2 shows the volumetric energy deposition in Pb using the 6 | 7 a 
activi ١ ۱ = 
boiling temperature, and is vaporized. This ablated material flows toward the sition and the energy densities of the Pb-Li after deposition are computed. | spectrum given in Fig. VI.4-1. The A*THERMAL (4) code was used in this calcu- o / 
| u . B ۳ 
rsity of 1 center of the cavity and intercepts the ions generated by the explosion. The This leads to a simple heat balance method of calculating the mass of vapor- lation. Note that 88 M or 22% of the target yield is in the form of X-rays. = p j 5 
energy associated with the ions is absorbed by the Li 7Pbg3 gas and does not ized Pb-Li. The amount of mass vaporized due to the X-ray energy deposition is com- "M. 0.1 2i LL. 
o Aim | 
directly impinge on the first surface. The gas heats up to a very high The energy deposition from X-rays can be described by exponential attenu- puted using a simple model. We calculate the total energy per unit volume x 1 1 | > 
— — i 
temperature and starts to release its energy by thermal radiation toward the ation inside the material. The photo-electric cross sections are those given required to vaporize the lead, 8.7 x 103 J/cm’. From the X-ray energy depo- z | J ۱ 0 
| l 1 l 1 
cool first surface. The liquid surface temperature increases upon receiving by Biggs. (1) Determination of the energy deposited from a given spectrum must sition calculation the depth of the material that will be vaporized can be o ol 上 bi p 
energy from this thermal radiation and condensation of the vapor. The higher also account for transport of any scattered photons. If, however, the primary determined. This adiabatic model becomes more accurate for very short X-ray k ۱ ۱ ۱ | Se S 7 8 
اس لب‎ = | 
surface t ature increases th d i j i i ۳ j j he ignored and an deposition times, i.e. instantaneous. : ind 0 
ur emperatu ase e vapor pressure and consequently increases interaction is the photo-electric effect, the scattering can g p 0-5 ۱۵-4 ۱-3 IO 2 ۱0-۱ | 0 IO 2 Lu 
the vaporization rate. The cavity pressure is very high after the initial X- exponential deposition profile can be assumed. The photo-electric effect com- We simply balance heat in the following manner. © 
PHOTON ENERGY (KEV) 
ray deposition, and varies by the combined effects of evaporation and conden- prises 90% of the interactions for photon energies up to 10 keV, 30 keV and 70 > 
— 3 3 
sation. It is of critical importance to beam transport and, to a lesser keV for carbon, iron and molybdenum, respectively. (2) Since lead is of yet Heat of vaporization of lead = 8.7 x 10° J/kg = 8.7 x 10> J/cm * 
3 ۱ 
Cp(Tb-To) Cp (1750°C - 400°C) = 1944 J/cm Fig. VI.4-1 X-Ray energy spectrum resulting from pellet burn. ul 
LU 


higher Z, the photo-electric effect is dominant to even higher photon energy. 







Total (J/cm?) for vaporization = 1.06 x 104 J/cm. 











degree, vacuum pumping to calculate the pressure history after an explosion. 










One important factor determining repetition rate is that the vacuum condition The volumetric energy deposition for a monoenergetic spectrum can be 
in the cavity be low enough to allow beam transport. represented by / 
In this section, the behavior of the Pb-Li vapor is discussed. In From the X-ray energy deposition curve (Fig. VI.4-2) this corresponds to a 
-u X -— | 
۰ thickness of approximately 2.5 microns. TO account for energy deposition 
F 






section VI.4.1 the deposition of target generated X-rays and the subsequent 8 v ov 
above Cp Tp - To) in the region between the saturated liquid and the sensible 




























































These calculations were repeated using cylindrical geometry and the re- 
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一 144 — 
M ۱ 
i ١ 
ul | heat, we have assumed that there is some additional vaporization. This addi- 
١ 
N | > tional vaporization accounts for the additional energy deposition above 
5 1 a or Cp(Tp - Tg) as is shown schematical!y in Fig. VI.4-3. This gives an addition- 
O 
N Q. | x al thickness of the vaporized layer of about 1.8 microns. Therefore the total 
d 
$ > | - thickness of material vaporized = 4.3 microns. 
١ 
oO 0 ١ — We now can calculate the mass of vaporized Pb-Li, 
م 8 لها‎ | a 
I | 
—* 8 م‎ ١ * 2 
* "m + م فى‎ M, = 6 
uJ ۱ d. 
لنا‎ | ۳1 | density of lead 
= 
- | 3 radius of the cavity 
| 
< | لیا‎ ô = thickness of material vaporized 
W ۱ = | | 
工 t Substituting, we get M, = 13 kg/shot. | 
0 
be 
< 


sults did not differ by very muc This is because the geometry is almost a 









square cylinder. 






An example of how the X-ray spectrum affects the mass of material vapor- 






ized is shown in Fig. VI.4-4. Three different blackbody spectrums were used. 


The softer 





In this example, the tota! energy contained in X-rays is 33.6 MJ. 
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The harder the spectrum, the more 





spectrum of 0.5 keV has a very short range. 










it penetrates through the material. We conclude that the mass of material 






vaporized will be determined by the X-ray spectrum and X-ray yield. This is 















determined by the specific target design. 
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VI.4.2 Pb-Li Gas Dynamics and Radiation 





As discussed in the preceding section, target generated X-rays vaporize 
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1g. VI.4-2. X-Ray volumetric energy deposition in liquid lead. 





several kilograms of Pb-Li off of the INPORT tubes. The resulting gas absorbs 










the target generated ions, reaches a temperature of about 1 eV and becomes 




















provide data tables of interna 
ty and radiation opacities of : 
temperature, gas density and r 
occur in either the SahaU) or 
follow a semi-classical formal 
density is calculated self-con 
considered with the 20 lowest 
states for each gas species. 
Some results of calculati 
Figs. VI.4-5 and VI.4-6. Figu 
Pb83Li17 over the plane of gas 
plot where the Saha model is u 
can be a function of density v 


between the two regions follow 


where kg Tgas is in ev. (8) Sé 
gas temperature. It is also ` 
ionization on the gas temperat 
ideal. 

Figure VI.4-6 shows the ! 
gas temperature - radiation ti 
features of this are the tran: 
and the extreme variability 0 


radiation transport very sens 







partially ionized. As this gas spreads from the tubes into the center of the 


reactor cavities, it radiates photons which put a heat flux back onto the 

















tubes. Eventually, the gas uniformly fills the cavity and cools enough that 





the heat flux on the tubes is insignificant. Before this happens, though, 
there may be some additional vaporization of Pb-Li due to the heat flux. As 
the gas radiates, it also cools, affecting the rate of recondensation. For 
these reasons, it is important to consider the behavior of this gas between 
the time of vaporization and recondensation. 


To study this problem, it is necessary to know the equation-of-state and 






















optical data for gaseous Pb-Li. We describe an atomic physics calculation 
providing this information in section VI.4.2.1. 


Once the atomic physics of the gas is understood, the motion and radi- 





ation of the gas can be studied. In section VI.4.2.2 we discuss a computer 


BLACKBODY SPECTRUM 





simulation of the gas behavior which uses the results of the atomic physics 





INTO FILM 


calculations in section VI.4.2.1. 


q(x) 
| 
| 
| 
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With this computer model, we can predict what heat fluxes the INPORT 
tubes might experience from gas radiation and we investigate the sensitivity 


of such radiation to the total mass of the Pb-Li vaporized. In section 


DISTANCE 


VI.4.2.2 we present heat fluxes and average gas temperatures versus time for 


Fig. VI.4-4 X-Ray energy deposition for different blackbody spectra. 


DISTANCE INTO LEAD FIRST WALL (microns) 


various Pb-Li masses. 
Finally, in section VI.4.2.3 we state our conclusion that the gas behavi- 


or is fairly sensitive to the mass of the gas and comment on the significance 


ENERGY DEPOSITION FOR 33.6 MJ X-RAYS 

















of these results. 


Fig. VI.4-3 Conceptual representation of Pb-Li vaporization. 


VI.4.2.1 Equation of State and Opacity of Gaseous Pb-Li 
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The atomic physics of gaseous Pb-Li is studied with the MIXER (5:6) com- 










puter code. This code has been developed at the University of Wisconsin to 
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VI.4-10 1 ١ 
د د فى فى ىك فى فى فا ف فى فا فى فا 3 لا ظ ظ‎ | 
| energy density, ionization state, heat capaci- VI.4.2.2 Radiation - Hydrodynamics in Pb-Li xo 8 P3 2 只 5 N 8 8 N 8 = = 8 > gi 
jas mixtures. This is done as functions of gas We have modeled the behavior of the Pb-Li gas with the 1-D Lagrangian | «e [TTTTTTTTTTTTTTTTTTTTTTTTTTTTTIT 
adiation temperature. Ionization is assumed to radiation-hydrodynamics computer code, FIRE. (10) This code has explicit | " 
the Coronal (8) model and radiative processes to hydrodynamics and implicit energy transport where heat is conducted through 7 E 1 
ism.(9) In a mixture of gases the electron two fluids - the gas at a local temperature of Tqas and the radiation fluid at "1 f 
sistently and the interaction of radiation is a temperature of Tradiation. In these calculations, conduction through the i = 
excited states of the six most common ionization radiation field dominates the thermal transport. Transport coefficients and , 
w E in ! equation-of-state information is obtained from a table of data provided by a er (0 
ons with this code for Pbg4Li,; are shown in E o 0 ۴ 7 calculation with MIXER. 1 
re VI.4-5 is a plot of the ionization state of ص‎ = $1 LE FIRE cannot accurately model the behavior of fluid which is as dense as a 1 
temperature and density. One can see on this 5 i N A liquid. Thus, we must analytically calculate the dynamics from the time when -r U 
sed by remembering the fact that Saha ionization he m 0 EE: the Pb-Li is at liquid density until the time that the density is low enough | h 
hile Coronal ionization cannot. The boundary * E A for FIRE calculations to be appropriate. We assume that the gas obeys an iso- | | - 
5 the criterion z s A 0-4 EE: thermal blow-off formalism where the temperature of the gas is that due to the | -[ « 
e = ` 222 deposited target generated X-ray energy minus the energy of vaporization and | J 
10° (kg — cm ١ the energy of ionization. We arbitrarily assume that the energy from target 1 - 
9 generated ions is uniformly deposited in the gas at 1.5 x 10-4 seconds after ۱ p” or 
ha ionization is used at high density and low vaporization and that FIRE can be used any time after this. Thus, at this [ 
important to note that the dependence of the time we start the computer calculations with the gas density profile being the i ١ 
ure makes the equation-of-state rather non- Gaussian shape predicted by the blow-off model and the total energy in the gas N p 5 
being the target generated X-ray and ion energy minus the vaporization and [ 
lanck averaged mean free path plotted over the ionization energy. ۰ b 
əmperature plane for a given density. Notable The FIRE calculation simulates the gas dynamics until the time the gas i p » 
sparency of cool gas to low temperature radiation reaches the center of the cavities. Figure VI.4-7 shows the positions of the [ 
f this parameter. This variability can make | PbassLil7 PbesLiiz s Lagrangian zone boundaries for a typical calculation during this period. Upon re q ۱ 
itive to the changes in the conditions in a gas. CHARGE STATE PLANCK MEAN FREE PATH reaching the center, the gas begins to convert its bulk kinetic energy into p اس‎ Maa ا‎ d 
| Fig. VI.4-5 a Pb DENSITY = 2.7 09 cm? heat, but FIRE, being a one-dimensional code, predicts that the gas will z © Y 5 2 D Y 品 号 e 2 e 8 o © 
: N N > + * n m N N N 一 مم‎ 
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Fig. VI.4-11 
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Fig. VI.4-10 
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Fig. VI.4-8 
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reflect off the center and propagate back towards the tubes. This is not 





physical because the system does not have the symmetry needed for such reflec- 





tions. For this reason we stop this first simulation calculation at the time 
when the gas reaches the center, typically about 1 millisecond after vapori- 
zation. 

To simulate the behavior past this time we assume that the gas has a uni- 
form temperature and density and has no bulk kinetic energy immediately after 


the gas collapses in the center. We find the temperature that will give the 


proper amount of energy, 


Egas ۳ x-ray 5 ion E -E 


Fvaporization ^ ‘ionization ^ ‘radiation ° 
where Eradiation ÎS the energy radiated back into the tubes by this time. We 
thien use FIRE again with this new initial condition and with the gas con- 
Strained not to move. By combining both uses of FIRE, we calculate the heat 
flux on the tubes due to radiation and the average temperature of the gas 
versus time. 


To test the sensitivity of the gas behavior to the total mass of Pb-Li 


vaporized, we have completed calculations for three different total masses; 





6.7 kg, 13.3 kg and 26.6 kg. In Figs. VI.4-8 through VI.4-10 the heat fluxes 





on the surfaces of the tubes are shown for total masses of 6.7 kg, 13.3 kg and 


26.6 kg, respectively. 





The average temperatures are shown in Figs. VI.4-11 through VI.4-13. 





Notice that there is in each case a large pulse of heat reaching the surface 










immediately after the gas collapses in the center. This is an effect of the 









conversion of bulk kinetic energy into heat. As is seen in Figs. VI.4-11 


through VI.4-13, there is also a pulse in the average temperature of the gas 
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Fig. VI.4-7 
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vI.5-2 n VI.4-22 — 160 — B ۷1.4-23 ee 
一 which is due to this conversion of kinetic energy into heat and which causes aT. oc, 31 (VI.4-1) The boundary conditions for 
= ey aes Pe eS pret = O the increase in heat flux. The same process has been seen in imploding wire 9X 
z 1 1 exper iments. Mond. ` 3.64 ۳ (M/RT 人 (VI.4-2) 
1 E ۱ | VI.4.2.3 Conclusions | 
5 N Comparison of Figs. VI.4-8 through VI.4-10 shows that the heat flux can Map. ^ P. rr (VI.4-3) K. 
ul P E 0 vary by two orders of magnitude when the total mass is changed by a factor of 
4 four。At 1ow gas mass the temperature must be high because the internal energy | a & is thermal conductivity of Li,7Pbg3 in which qtot = qr + Qc = 0 
< ۱ 5 density is high. Also, the radiation passes through less dense gas more easi- : 9 is density, g/cm? 
" f |o ly. Thus, the heat fluxes are higher for low total gas mass. As the mass of 3 is specific heat, J/g K qp is energy deposited 
~ P 3 - 1 0 3 the gas is increased, both of these trends are reversed and the gas radiates 1 is temperature, K qç is energy deposited 
- A o : 4 | energy more slowly back to the tubes. ۱ اه‎ fran ticot surfati GÀ رو‎ is energy loss due i 
1 2 N - 3 t- VI.4.3 Evaporation/Condensat ion Calculations : t is time, sec h is heat transfer coe 
> > O ۲ 1 t $ Once the radiation heat flux and the temperature of the gas is known, the — is condensation rate, g/sec - cm T, is temperature of ti 
= 9 2 = 1 0 E behavior of the gas density in the cavity over time can be investigated. The | Pa is pressure of the gas, Torr ظ‎ 
wo ١ : rate that gas leaves the cavities is governed by geometry of the cavities and M is molecular weight | T(x) is the temperatur 
E 1 1 the gas temperature。 The rate that gas is added to the cavity through evapor- | J is temperature of the gas, K ture is unknown at the begii 
Wa 0 b o ation off of the tubes is determined by the temperature of the Pb-Li film on nevp. is evaporation rate, g/sec - emê | However, after a few cycles 
O P 1 9 the tubes。 This temperature can be found through a temperature diffusion P. is saturation vapor pressure of first surface T(y) can be calculated. 
N $ - calculation in this film for the surface heat fluxes calculated in section 1 is the temperature of the first surface, ۰ Equation VI.4-1 cannot 
4 ۱ ۱ VI.4.2. | and boundary conditions. T 
- ۲ ۱ In section VI.4.3.1 we describe the methods used in the temperature The initial conditions are: used: 
= diffusion calculation, in the condensation and in the evaporation. In section 
VI.4.3.2 we present and discuss results. Conclusions of a general nature are p T= 11501 0 > x < 1074 cm T; 5 
١ made in section VI.4.3.2. T = T(x) RT 1074 cm 
(As 2031 5170 393۷ VI.4.3.1 Calculation Procedure d 7 - 108 Torr 
The governing equations for heat transfer, condensation and evaporation Ty = 14500 K in which A6 is time step an 
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Fig. VI.4-22 Temperature of the first surface versus time for different 
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Fig. VI.4-16 Condensation and evaporation rates on first surfi 
of vaporized gas = 13000 g. 
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cycles. 


FROM THE FIRST 
(cm) 


SURFACE. 


DISTANCE 


j 
0.075 0.105 0.135 


L 
0045 


L 





۱ 
O 5 


300 ۱ 


— 179 一 


Bap "Uv 






CYCLE 3 


CYCLE 2 





| 
0 
0 
t 


500 


(De) 1178-1 


VI.4-41 


Mo = 13,0009 











VI.4-26 0 
N 
o 
4 = 
f 
~ Q 
* 
42 
t 
ل‎ 9 
u 
ده‎ 
£ 
© w 
= > 4 2 
3 
42 
42 
۳ ۳ 
۹ 9 ۴4 
ها‎ 
له‎ 
<= 
> E 
j 9 c 
P c 
LJ 
به‎ e 9 
Q- 6 - هھ‎ * 
* 0 
72s © 
= — 58 
" wu eco 
LI g 
' +? 
d © w it 
- £ 
- N 
ها‎ o 
2» 
9 2 O 
4 ته و‎ 
هه‎ ce 
O & 
o & 
: $ 
co حب‎ 
= © 
w 
f 
a 
e 22 
و‎ 
PP 
3 + 
4 = as 
2 
et 
~œ سم‎ 
ل‎ ‘o دا‎ 
© 
4 1 1 1 1 E 4 e 
E BÀ BÀ D 
© H © ص‎ 
373144 
— 178 — 
VI.4-40 


(Se) 





M, = 13,000 g 








1100 


END ww 


before time step A6, cond; is the conductance e between zones j and i, 
[vec J; is the thermal inertia of zone i and 1 is new temperature of zone i 
after A6. 

By using Eqs. VI.4-2, VI.4-3, and VI.4-4, and starting from the initial 
condition, all the information in the cavity can be calculated step by step. 
VI.4.3.2 Results and Conclusions 

We have calculated heat transfer in the Pb-Li film, evaporation and 
condensation with the following assumptions: 

a. The cavity is spherical with R = 5 m. 

b. The gas zone has a uniform temperature and pressure. 

c. The condensation and evaporation are by Pb only. 

d. The material properties are independent of temperature. 


e. The first surface forms a spherical enclosure over the cavity. 





f. The thickness of the tube wall is 1.5 mm. 

g. The coolant temperature is 324?C and has a heat transfer coefficient 

of 12 W/cm*-*C. 

The results of a series of calculations for 13 kg of vaporized Pb-Li are 
summarized from Fig. VI.4-i4 to Fig. VI.4-20. Figure VI.4-14 shows the 
temperature of the first surface as a function of time. An examination of 
this curve shows that the temperature history is not the same for the first 
two shots. Figure VI.4-15 shows the temperature of the wall as a function of 
distance from first surface just before the first three explosions. It can be 
seen that temperature has almost converged to a quasi-steady-state after 3 
shots. Figure VI.4-16 shows the condensation and evaporation rates. Note 
that for a vaporized mass of 13 kg, evaporation is not important. Figure 
VI.4-17 shows the total condensation and evaporation and the mass in the 


cavity as a function of time. At times after 0.1 sec, the mass of gas in the 
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Eq. VI.4-1 are: 


dT _ 
=hT-T) at x= x, 


)。 Me have used the following definitions: 
due to radiation 
due to condensation 
,0 evaporation 
Fficient 


e coolant. 


ə profile in the coolant tube wall. This tempera- 
nning of the calculation and has to be estimated. 


, a quasi-steady-state will be established and 


be solved exactly due to the complicated initial 


nerefore, the following finite difference form is 


cond. . 


b ۴ 7 VI.4-4 


[voc, J. 
1 15 > TU T; is temperature of zone 1 
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Fig. VI.4-17 Total condensation, evaporation and mass in cavit 
mass vaporized = 13000 g. 
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first surface and, consequently, causes a large evaporation rate. Thus, the 
effect of evaporation dominates cavity mass transfer. Figure VI.4-22 shows 
the surface temperature of the three cases. As expected, the surface tempera- 
ture increases drastically if the mass of vaporized gas is reduced. This 
occurs because the heat flux shown in Fig. VI.4-21 is so much higher at low 
mass. This leads to a higher evaporation rate in the low mass (6.7 kg) case 
shown in Fig. VI.4-23. Here the evaporation rate is shown to be, at times, 
much larger than the condensation rate. This leads to the increases in cavity 
pressure, mass and atomic density with time at low mass shown in Figs. VI.4-24, 
۷1۰4-25 and VI.4-26, respectively. Fig. VI.4-27 shows the atomic density at 
0.2 second for the three cases of initial mass. There is a sharp minimum at 
Mo = 13 kg which is due to the gas temperature being high enough for quick 


condensation and the heat flux being low enough for insignificant evaporation. 


When M, - 26.6 kg the gas temperature is too low and recondensation is slow. 
When M, = 7 kg, the radiation heat flux is high enough that evaporation is a 
problem. 


In summary, we have considered the vaporization, gas dynamics and conden- 
sation of Pb-Li from the film on the INPORT tubes. The mass of the gas vapor- 
ized by the target generated X-rays is 13.3 kg and the gas density in the 
cavity 0.2 seconds after the target explosion is low enough for ion beam 
propagation. The gas density at 0.2 seconds has been found to be sensitive to 


the mass of vaporized Pb-Li and thus, sensitive to changes in the target X-ray 


spectrum. 





ities is a small number which is the difference of two larger numbers, M, 

| the condensed mass. Figure VI.4-18 shows the surface heating rate due to 
liation, condensation and evaporation, the summation of which is the heating 
the first surface. Notice that the heat from radiation is dominant for 

-h of the time. Figure VI.4-19 shows the pressure and Fig. VI.4-20 shows 
"ticle density in the cavity as a function of time. 

The vacuum pressure at time - 2 second should be low enough to allow 
am propagation. The particle density at this time is 4 x 1010 cm-3。 The 
lculation of the beam stripping cross section will show that this pressure 

acceptable. It must be noted, however, that the total mass in the cavity 
only 7 x 107? g. This is obtained by subtracting the condensation mass 
om the sum of initial mass and evaporation mass, both of which are large and 
ry close. The total condensation mass is 1.5 x 104 gm. Therefore, the 
certainty of the total mass in the cavity at t = 0.2 second is very large. 
change of assumptions may lead to a significant variation on the final 
rticle density. 

A change of the target design may change the energy associated with the 
.rays and/or the spectrum of the X-rays. This will change the initial mass 
iaporated by the X-ray deposition. To study the effect of different X-ray 
۱6۲۵۷ or spectra, a different initial mass is assumed. Two additional com- 
iter calculations are performed with M, - 26 kg and 6.5 kg. The larger mass 
xduces the initial gas temperature and also the radiation heat flux to the 
irst surface. Consequently, the thermal energy deposited on the first sur- 
ace is spread over a longer period of time, as can be seen on Fig. VI.4-21. 
he lower gas temperature also reduces the condensation rate, as can be seen 
n Eq. VI.4-2. Smaller initial mass, on the other hand, increases the initial 


as temperature. The high gas temperature results in a large heat flux to the 
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Fig. VI.4-22 Temperature of the first surface versus time for different V1.4-37 — 175 — 
initial cavity gas masses. 
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Fig. VI.4-27 Final atomic density in the cavity versus mass of PbLi 
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Fig. VI.4-26 Atomic density versus time for different masses of vaporized gas. 
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Fig. VI.4-25 Total mass in cavity versus time 
vaporized gas. 
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24 Cavity pressure versus time for different masses of vaporized gas. 
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rectangular ducts of the same dimensions, leading to cryopumps located in the 















































VI.5 Vacuum System 
-ember | boundary layer develops on the cavity surfaces. For this reason the gas 
| 1.5.1 General Description behaves more as in molecular flow rather than viscous flow. The behaviour of upper corners of the cavity as shown in Fig. VI.5-2. These pumps are well out i 
The previous section described the behaviour of the PbLi vapor in the : 3 i f 
od | | d ۱ the non-condensable species within the vacuum ducts will thus be treated with : of the radiation field and Monte Carlo analysis has shown that the effect of 1 
۲ ۲ ۱ cavity after each shot. It was found that the vapor pressure of PbLi falls | i iud 
۰ molecular flow theory rather than viscous flow. ' 全 neutron streaming on them is negligible. 
below 107^ torr in the order of 200 ms after the shot. In the present section 3 ۱ ۱ 
: Most of the recent fusion reactor design studies have utilized compound 2 VI.5.2 Pumping of Non-Condensable Species 
we will deal with the non-condensable fraction of the gas load within the e Í 
| | 1 i cryopumps for pumping hydrogen and helium. In these pumps, hydrogen species 3 2 Table VI.5-1 gives the source and species of the non-condensable gas 
cavity, namely Do, To and He. We will attempt to access the effect of the "m ۱ aa ۱ 1 N 5 5 
are condensed on chevrons coo 0 ~ 116 the helium is cryotrapped on oc e 
PbLi vapor on the vacuum systems for the cavity and the beam lines. * | 7 a MP i ۱ £ ~ 10.000 K. the 
ne | | molecular sieves applied to panels which are at 4.2 K. Typically these pumps á o ; At the initial gas temperature in the cavity (1 eV) of ~ 10, , 
itions The cavity atmosphere is primarily dictated by the requirements needed to ; 2 سا‎ i 5 
1 have speeds of ~ 5 t/s cm" for Do, To and ~ 2 ۸/5 cm for He. A major 015- 2 non-condensables constitute a gas load of: 
propagate the beams to the target with minimal losses due to stripping and = 
Aw | advantage of any cryopump is the need for periodic regeneration. This is done t 
- charge exchange. These issues have been addressed in Chapter V. Similarly, 7 d 5 
| | by valving the pump off and warming up the cryosurfaces. Although regener- 3 D = 320 torr. £/shot 7 
| the equilibrium pressure of the non-condensable fraction has to be maintained A bet Í ul 
- | ation may only take - 20-30 min., during this time the pump is not available o ۲2 = 265 torr. £/sho z 
reasonably low, although their effect on beam propagation is not as great as a | | 5 
| T : | for operation. Some designs have resorted to providing twice as many pumps as 8 He = 326 torr. 2/shot | ui 
PbLi at the same number density, simply because of the difference in mass. | , v b» tê 5 ۰ ۱ 5 ۰ 2 | | z 
۱ ۱ ۱ needed, such that only half are on line at any one time. s is not space T | " 
h Furthermore, the noncondensable partial pressure has to be kept low because it 1 l ۱ 
و‎ A | p | , . P ۱ efficient, particularly where there are space limitations. In several past 5 o The vacuum ducts are assumed to be about 400?C in order to condense any PbLi N 
constitutes a continuous source of molecules migrating into the beam lines * ۱ ۱ ۱ 
2i 7 i ۱ : designs(1,2) the UW group has proposed the use of back to back cryopumps such à n vapor that they admit. At this temperature the throughput becomes: 
-— where the pressure must be kept at ~ 10 torr. | ox مخ‎ ۹ i. 
| of | | | as that shown in Fig. VI.5-1. Two sets of cryopanels, with integral baffles SS kh = p 
小 | The source of hydrogen in the cavity is from the unburned fuel in the | Fa g 2 2 e 
| | and chevron shields, are provided back to back. While the front panel is in z 三 à - - D> = 108 torr. £/sec. 
target, the newly bred T released within the cavity and the deuterium gas load e 
use, the rear panel is being regenerated and the pump body itself constitutes 3 3 < S ۲2 = 89 torr. £/sec. 
from the target injector. Helium comes primarily from the DT reaction and 8 oo » 
the shut-off valve. Although some development work will be needed to perfect 5 P — He = 110 torr. £/sec. 5 
To breeding, with a minor contribution from To decay. At pressures on the : 5 ۳ ۱ ' نسح‎ to buio 00 leo for ediutenttelty «j^ E ۳ uj 
e seals for such a pump, we feei that it holds the promise for su "t | ۳ 
order of 1074 torr, outgassing does not dominate and will be neglected.  Out- Win ^ A Y Y ۱ ۱ ۱ lar fl is given by the ó 
۱ " increasing the pumping capacity in systems which are space limited. We are n ^ ^ ۱ The conductance of a short duct ) < 25 r) for molecular flow 15 9 y " 
gassing is a serious problem in the beam lines where the pressure must be e ^ A Y Q ١ الو‎ : De / 
| proposing such pumps for the present study. | is N A Y ) | ] ^ relationship: «z 
three orders of magnitude lower. | ۱ i S A 2 Y B $ 
‘ In the HIBALL design, vacuum ports have been provided at the junction m ^ A ۱ 7 2 5 
Immediately after a shot, the walls of the INPORT tubes as well as other | o SY K | | / 1 3 LU, 1 ir sec — 
| between the upper blanket modules and the top of the INPORT tubes. The spaces 4 سم / ول ال‎ (TE =y + x) RT —3 
cavity surfaces act as condensers for the PbLi vapors. Although the pressure duct A cm 
between the headers which connect the two blanket systems together are devoted g ۳ 
reaches ~ 100 torr, because of the sticking coefficient of unity for PbLi, no g 
1 to vacuum ports. There are 30 ports 65 cm high and 120 cm wide, connected to su 
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